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Abstract 

Some quantum properties of QED 3 are studied with the help of an exact evolu- 
tion equation of the effective action with the bare fermion mass. The resulting 
effective theory and the occurrence of a dynamical mass are discussed in the 
framework of the gradient expansion. 



1 Introduction 

QED in 2+1 dimensions (QED3) was introduced as a toy model for chiral 
symmetry breaking [1, 2]. There the Schwinger-Dyson equation for the fermion 
propagator was solved, in order to study the dynamical fermion mass gener- 
ation. More recently, similar studies were done analytically and numerically 
[3] , where different resummations are considered in the study of the Schwinger- 
Dyson equation, taking into account different ansatze for the the vertex. 

Besides the toy model motivation, QED3 happens to be relevant for the 
effective description of planar antiferromagnetic systems [4] . Indeed, the d— wave 
pairing that occurs in these systems allows a linearization of the dispersion 
relation around the nodes in the energy gap, which leads to relativistic-like 
effective theories. 
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Motivated by these recent developments, an alternative approach to the 
quantum properties of QED3 is proposed here. This approach is based on an 
exact functional method, where the quantum fluctuations are controlled via the 
bare fermion mass. When this mass is large, the interaction can be neglected 
and the quantum corrections are weak: the system is almost classical. As the 
bare mass decreases, the quantum fluctuations gradually appear and the system 
sees the generation of quantum effects, leading to the the full quantum effective 
theory. 

This functional method was also used in 3+1 dimensions for a scalar theory 
[5] and for QED [6], in the description of a coulomb gas [7] and in a Wess- 
Zumino model in 2+1 dimensions [8]. In all these cases, the usual renormalized 
parameters were recovered at one-loop. Beyond one-loop, the resummation 
provided by the gradient expansion used in this functional method docs not 
coincide anymore with the loop expansion. In the example of QED and for the 
one-loop approximation, a singularity appears in the coupling constant flow as 
the bare fermion mass decreases, i.e. as the fluctuations appear. This singularity 
is the analogue of the famous Landau pole and it disappears beyond the one- 
loop approximation, i.e. when taking into account the full equations obtained 
in this method. 

In section 2, we set up the frame of this method and give the resulting flows 
that are obtained, the computations being given in Appendix A and Appendix 
B. We show then in section 3 that the resulting effective couplings are consistent, 
at one-loop, with the usual renormalized parameters computed with standard 
Feynman graphs. In section 4, we give the non-perturbative solution and discuss 
the generation of a dynamical mass in the massless case ("massless" refers here 
to the bare theory). A self-consistent relation is found between the effective 
coupling and the dynamically generated mass. This relation shows that the only 
way to have a vanishing dynamical mass is to have a non-interacting effective 
theory. This result is set in the lowest order of the gradient expansion and thus 
the influence of the number of flavours is not discussed, as explained in section 
5 which is let for a discussion. 

2 Evolution equations 

We consider a parity conserving theory, what is possible with 27V fermions, 
constituing N couples in which the masses are opposite. The 2N fermions are 
then written as N 4-component fermions and we consider the 4x4 reducible 
representation of the 7- matrices given in [2] . 
The starting point is the bare action 
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where T^ v and L^ v are respectively the transverse and longitudinal parts of the 
photon propagator and a is a usual gauge-fixing parameter. 

The control of the quantum effects consists in playing with A. As this pa- 
rameter goes to oo, the fermion bare mass is large and thus freezes the quantum 
fluctuations: the system is classical and the effective theory is the bare one. On 
the other hand, as A decreases, the quantum corrections gradually appear in the 
effective theory. 

The idea of the method is to compute the evolution of the effective action T 
(the proper graphs generator functional) with A. Starting from A = oo, where 
r = S, the massive case will be obtained when A decreases down to 1 (see 
section 3) and the massless case when A decreases down to (see section 4). 

The evolution equation for T was derived for N — f in [6] and we review in 
the Appendix A this derivation for any N. The result reads 



N 
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d x r = m B Tr\ i-^k + i > ( 2 ) 



and we recall that this evolution equation is exact and does not rely on any 
approximation. In the framework of the gradient expansion, we then make the 
projection of T on the following functional subspace: 



T = J d 3 x{^((3 T T^ + af3 L L^)A v 

N 



(3) 



where (3t, Pl, z, g,m are functions of A. The step (3) constitutes our approxi- 
mation and is valid in the IR. It is shown in the Appendix B that plugging this 
ansatz in the evolution equation (2) leads to the following differential equations: 



dxfi L = (4) 
d\ m = tub 

OxPt = -N- 
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It is interesting to note the following points: 

• As expected, the photon does not acquire a mass in the evolution with A. 

• The longitudinal part (3l of the photon propagator does not evolve with 
A. 

• In the Landau gauge (a = oo), there is no fermion wave function renor- 
malization and z = 1. 

• The fermion mass m has a classical flow with A (no quantum contribution). 
As a consequence, we have 

Pl = 1, (5) 
m = ms(A + c), 

where c is a constant of integration linked to the dynamical mass generated in 
the massless case since. 

m dyn = lim (— ) = c— , (6) 

where z = 2(0). 

If we note @t = (3, the remaining differential equations can be written 
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Finally, we will consider the Landau gauge where z = 1 and the differential 
equations that remain read then 

B a ^ 1 

Ay 97r/3m B (A + c) 2 ' 

In this case, the dynamical mass is given by 

mdyn = cm B . (9) 
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We remind that the bare fermion mass is Xms and not tub, such that the 
massless case is reached when A — ► 0. tub must be seen as an initial parameter 
which cannot vanish (the whole method is based on ^ 0!). The possibility 
of having md yn = then occurs if c = 0. 



3 One-loop solution (massive case) 

To obtain the renormalized parameters in the massive case, we integrate the 
differential equations (8) from A = oo where the quantum corrections are frozen 
(3(00) = <7b,/3(oo) = 1), to A = 1 where all the quantum effects are generated 
(we note g(l) = fli,/3(l) = ft). 

Before going to the full treatment of the differential equations (8), let us 
check that these lead to the expected results at one-loop. We note for this that 
the right-hand sides of Eqs.(8) are proportional to h, such that the one-loop 
approximation consists in replacing the parameters of the right-hand sides by 
the bare ones (and thus take c — > since this constant arises from quantum 
effects). We have then, restoring the factors U, 



™ - ~ nN Sr B h + - (10) 



QirrriB A 



where the dots stand for higher orders in K. The integration of these equations 
from A = 00 to A = 1 gives then 



ft - l + %N-^- + ... (11) 
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It can be checked with standard computations that these results are given by 
the one-loop Fcynman graphs. Indeed, the graph corresponding to the one-loop 
vacuum polarization (for one fermion flavour) is 
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where q — p — (1 — x)k and qs denotes the Euclidean momentum. As expected, 
the would-be linearly divergent integral (second line) cancels when using a gauge 
invariant rcgularization. The one-loop function (3\ for TV flavours is then given 
by Eq.(ll). Note that this result is consistent with the limit k « tub of the 
large- N polarization tensor [1]. 

The graph corresponding to the one-loop correction to the coupling is, in 
the Landau gauge and in the limit of zero incoming momenta, 

J p V V -m B P m B 

JA^p f 1 Pl + 3m| 
9 7 J PE p 2 E (p E + mlr 

= — 

which leads to the result in Eq.(ll). 

4 Non-perturbative solution 

We turn now to the complete resolution of the differential equations (8). Taking 
the ratio of these two equations and integrating, we have for any A 

P=( 9 A • d4) 
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Plugging this result in the differential equation (8) for g gives 

-3JV/2 1 

The integration from A = oo down to any A leads to 

9{ X )=9B 1 + . (16) 



18tt m(A) 

The solution for the photon propagator renormalization is then 

= ( JUL V"' 2 = (i + 2^±ij* N 3 "" 3A ' +4) 



={ 1+ —^>) ,17) 

So as to discuss the dynamical mass generation, we should take the limit 
A — > 0, where m(0) = md yn and w e note g(0) = g - We have then 
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From this last equation we see that the only possibility to have nidyn = is 
go = 0. It is indeed expected that no dynamical mass can be generated in 
a non-interacting effective theory. Eq.(18) should be seen as a self-consistent 
relation between the renormalizcd parameters go and md yn , relation which is 
parametrized by g B . This result, valid in the present approximation of the 
gradient expansion, is independent of N. 

5 Discussion 

Let us first show the link between the equations obtained by this method and 
the Schwinger-Dyson equation for the fermion propagator. 

One can see from the Appendix B that the evolution equation of the fermion 
propagator is 

dxG-\k) = -m B -im B g 2 z 2 I lli G 2 {p) lv D^{k+p) 

Jp 

= -m B -im B g B / l^G 2 B {p)-i v D B v (k+p) 
Jp 

+ higher orders in h, (19) 

where G B and D B V are the bare fermion and photon propagators respectively. 
But since G B (k) = tf — \m B , one can note that 

d x G B (k) =m B G 2 B {k), (20) 

such that the one-loop approximation to the evolution equation of G _1 can be 
written 

d x G-\k) = dxG B \k) - ig% f ^d x G B ( P h„D B u (k+p). (21) 

Jp 

Since D B V is independent of A, the previous equation integrates as 

G-\k) = G B \k) - ig% f lti G B {p) lv D B »{k+p). (22) 
Jp 

Therefore the one-loop approximation to the evolution equation of the fermion 
propagator gives the usual one-loop self energy, which is also the one-loop trun- 
cation of the Schwinger-Dyson equation for the fermion propagator. The higher- 
order terms in % are different though, since the resummations provided by the 
two methods are not the same. 
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The study of dynamical mass generation with the Schwinger-Dyson method 
[1, 2, 3] takes into account the momentum dependence of the fermion and photon 
self energies, what is not the case with the gradient expansion truncated to the 
present IR approximation. For this reason, the influence of the number of 
flavours is not discussed here. Indeed, it is the momentum-dependence of the 
fermion self energy which leads to a condition on the number of flavours: it was 
found in [2] that there is no dynamical mass generated above N c = 32 /ir 2 , and 
the studies made in [3] agree with the existence of N c , and of the same order. 

The next step thus, to discuss the influence of the number of flavours, would 
be to go further in the gradient expansion, so as to take into account the mo- 
mentum dependence of the self energies. But this would lead to "partial-integro- 
diffcrcntial" equations (derivative with respect to A and integration with respect 
to the momentum) and would involve new approximations, similar to the ones 
used in the study of Schwinger-Dyson equations. The present scheme is based 
on the gradient expansion approximation only and is therefore believed to give 
a relevant insight into this problem of dynamical mass generation, besides the 
study of the role of N. 

To obtain an explicit expression of the dynamical mass as a function of the 
bare coupling only (and therefore compute the value of the constant c in Eq.(9)), 
one would need to take a finite initial value for A, corresponding to a system 
where some quantum fluctuations are already generated. The initial values of 
the parameters defining the theory could then be taken from a loop-expansion, 
valid when A >> 1. But the result would then depend on the initial value of A, 
whereas the motivation of the present work was to find a self-consistent equation 
like (18), independent of the A-scheme used here. 



Acknowledgements: I would like to thank Janos Polonyi for useful discus- 
sions. 

Appendix A: Evolution equation for V 

The bare action, functional of the fields ip, x/j, is 



S = J d 3 x{^A»(T^ + aL^)A» (23) 

JV 

k=i 

and the connected graphs generator functional, function of the sources ?7fe,?7;,i M 
is defined by W = — i\nZ where 
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Z = J P^vUjexpjiS + i J d^^i^ + ^^fe+^fe^j. (24) 
The functional derivatives of W define the expectation value fields ip k , tpi, A^: 



SW 1 ? 

0V1 2 



5W _1 - 

$r] k Z 
8W 1 
~Z 



-<A^>= Ap, (25) 



where for any 



< 4> >= J V[4>] exp |*5 + d 3 .x + J^^kVk+rjJk^j j . (26) 



We also have 

Inverting the relations (25) which give the expectation value fields as functions 
of the sources, we define the Legendre transform r of W by 



T = W- J d 3 x^A^ + J2^ k Vk+vM'j. (28) 
From this definition we extract the following functional derivatives: 



5T_ 

Sipi 
6T 

ST 
5A^ 

S 2 T ( 5 2 W 



= Vi (29) 
= -Vk (30) 
= ~U (31) 



The evolution of W with the parameter A is, according to (27), 
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dxW = ~ f d 3 xf2<ki>k> 
J fe=i 

/N . N £2yy 
d 3 x il) k il) k - im B j ^X^^F ' ( - 32 -' 
/ — i J i — i VkVk 



k=l ° fe=l 



To compute the evolution of T with A, we have to keep in mind that its inde- 
pendent variables are tp k , tp^A^ and A, such that 



f , \ SW v\ / SW SW 

- j d 3 x IdxfA^ + Y, @kdxVk + dxrj k ^k) j 



fe=i 

= d x W. (33) 

Combining these different results, we finally obtain the exact evolution equation 
for the proper graphs generator functional T: 

W-*//.t|-i* + (( j g ; )" 1 ) (34) 
= m B 5 3 (p + q)Y\-' l l ) k{p)^k{q)+i 
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where we denote 



5 2 T 



Hkip)Hk(q) 



/(..., = /^(...,. 

Appendix B: Evolution equations for the param- 
eters 

We decompose the second derivative of T in a sum of the field-independent 
diagonal part A and the part which contains the fields. The inverse of S 2 T is 
then taken as the following expansion 

(^r)" 1 = A -1 - A _1 ©A _1 + A _1 0A _1 6A _1 

-A- 1 0A- 1 eA- 1 eA- 1 + ... ) (36) 
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and this expansion is truncated up to the desired order consistent with the 
approximation (3) that we consider, i.e. up to the order 3 in the non-diagonal 
part 0. The result is then (the flavor indices are not written) 



-29 ll D^9 u G(p-<l>G(f>G(f> + ...)g}, (37) 

where (in Fourier components) 



D^(P,Q) = -^(^rT^(p) + -]-L^(p))5 3 (p + q) 



p 2 \(3 T a/3. 



L 



<f>(p,q) = -zg 4(-p-q) 

Pfaq) = -zg^(-p-q) 

f{p,q) = zg^(-p-q)Y, (38) 



with 



L,M = ML. (39) 
The contribution of Eq.(37) to the fermion propagator is then 

Tr \gTd^6 u g} (40) 

= -g 2 z 2 Jj}(k) [J rG 2 ( P )^D^(k + p^j ip(-k) 

^z 1 1 f- 

^>(fc) tfi/)(—k) + higher orders in k, 



6irm 2 a@L Jk 

where the higher orders in k are not considered in the present approximation 
(3) of the gradient expansion. We note that the mass term (order k°) vanishes. 
The contribution of Eq.(37) to the photon propagator is 



Tr {G(j)G(j)G} (41) 
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/ A„(-fc) ftr f G 2 { P )rG{k+p)AA v {k) 



[ A^{-k)A v {k)k 2 T^{k) + higher orders in k. 



As expected, no photon mass is generated and the correction to the photon 
propagator is transverse. Finally, the contribution of Eq.(37) to the vertex is 



+ higher orders in k and p. 

The last contribution G<pG4>G<f)G in Eq.(37) gives a cubic term in the photon 
field and vanishes by symmetry. The identification of the right-hand side and left 
hand side of the evolution equation (2) leads then to the differential equations 
(4). 

References 

[1] R.Pisarski, Phys.Rev.D29: 2423 (1984). 

[2] T.Appelquist, M.Bowick, D.Karabali, L.Wijewardhana, Phys.Rev.D33: 
3704 (1986). 

[3] P.Maris, Phys.Rev.D54: 4049 (1996); V.P.Gusynin, M.Reenders, 
Phys.Rev.D68: 025017 (2003); See references therein. 

[4] R.Shankar, Phys.Rev.Lett.63: 203 (1989); N.Dorcy, N.Mavromatos, 
Nucl.Phys.B386: 614 (1992); K.Farakos, N.Mavromatos, Phys.Rev.B57: 
317 (1998); M.Franz, Z.Tesanovic, O.Vafek, Phys.RevB66: 054535 (2002); 
I.Herbut, Phys.Rev.B66: 094504 (2002); See also references therein. 

[5] J. Alexandre, J.Polonyi, Ann.Phys.288: 37 (2001). 

[6] J.Alexandre, J.Polonyi, K. Sailer, Phys.Lctt.B531: 316 (2002). 

[7] S.Corrcia, J.Polonyi, J.Richert, Ann.Phys.296: 214 (2002). 

[8] J.Alexandre, hep-th/0306039, to appear in Phys.Rev.D68. 




12 



